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Abstract. We study structure functions of rotation measures in the Canadian Galactic Plane Survey (CGPS) region 
and the North Galactic Pole (NGP) to extract the interstellar medium (ISM) fluctuation information. The CGPS data 
are divided into three longitude intervals: 82° < « < 96° (CGPSl), 115° < I < 130° (CGPS2) and 130° < I < 146° 
(CGPS3). The structure functions of all three regions have large uncertainties when the angular separation is smaller 
than (50 fa 1°. A power law can fit the structure function well for S6 > 1°. The power law indices get smaller from 
CGPSl to CGPS3 and the amplitudes decrease. The variations of the large-scale field and the electron density have 
only negligible effects on the structure function and thus cannot account for the changes, indicating that the turbulent 
properties of the Galactic ISM are intrinsically longitude-dependent. The Kolmogorov-like fluctuation spectrum of the 
electron density or the magnetic field should produce a power law structure function with an index of 5/3 or 2/3, neither 
of which is consistent with our results of small indices in the three sub-CGPS regions. For the NGP region, the structure 
function is flat, showing that the rotation measures are mostly intrinsic to the extragalactic sources, and the ISM is 
very random in that part of our Galaxy. It is obvious that the ISM fluctuation is latitude-dependent when comparing 
the results in the NGP region and the CGPS regions. 

1 Introduction 

The electron density irregularities have been studied either by refractive scintillation or by diffractive 
scintillation phenomena as systematically summarized by Armstrong et al. (1995). The fluctuation 
spectrum from scales of ^ 10^ m to scales of ^ 10^^ m can be fitted by a single power law with a 
power index of ~ 3.7, i.e. the Kolmogorov spectrum. This spectrum can be extended up to 10^^ m. 
The property of magnetic field fluctuations is not yet clear although they are more important than 
the electron density fluctuations in some circumstances such as in the solar wind (Armstrong et al. 
1995). 

The structure function of rotation measures (RMs) contains information on electron density and 
magnetic field fluctuations (Simonetti et al. 1984). More data are available now allowing a better 
determination of the structure function and a careful comparison of the turbulent properties of different 
regions. In this contribution, we show the RM structure functions of extragalactic sources in the 
Canadian Galactic Plane Survey (Brown et al. 2003, CGPS hereafter) and briefly discuss the results. 
The detailed report of this work will be given elsewhere. 



2 The Theoretical Structure Function 

The RM structure function (-Drm), if assumed as a stationary random process, can be written as 

DnuiSe) - {[RM{eo) - RM(0o + S9)]^) , (1) 

where 69 is the angular separation of two sources for RMs in degree throughout this paper, and (■ • ■) 
denotes the ensemble average over any 6*0 so that the result is independent of Oq. Both the electron 
density (ng) and the magnetic field (B) are assumed to consist of uniform backgrounds {uq and Bq) 
and Gaussian fluctuations (Sn and JB) with averages of zero, in the form of 

rig = uq + 5n , 
B = Bo + (5B . 
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The exact spectra of the electron density and magnetic field fluctuations are not known, so it is 
straightforward to take the power law form as 

where q is the wave number, and C% are the fluctuation intensities. Ion and IqB are the outer 
scales, and a and (3 are the power indices. The structure function can be reduced to the form below 
when Ion = loB = lo = ^/Qo (Minter 1995) 

Dnum = ^RM{S9)+2Cl^BlClf{a)L"-'69"-^ + 2Cl^nlClg{p)L''-'59''-^ 

+ 2Cl^ClClh{a, /3)L(«+/3)/2-lg3-(a+/3)/2^^(a+/3)/2-2 ^ 

where ^BMiSO) is the geometrical term attributed to the variation of the large-scale electron density 
and the magnetic fleld, L is the path length, f{a), g{l3) and h{a,j3) are constants containing F- 
functions. As 5^ ^ 1 rad, the geometrical term can be approximated as 

*rm(^6I) « Cl^nlBlL^Se^ « 3 x IQ''^ {RMf 60^ , 

here 59 is in degree. If both the electron density and the magnetic field fluctuations follow a Kolmogorov 
spectrum, the structure function has a simple form as below, 

where C is a constant related to the intensities in the regular and fluctuating components. 
3 Analysis of Data 

In principle, two methods can be employed to extract the fluctiiation information of the ISM from 
RMs: the autocorrelation function and the structure function. The autocorrelation function is the 
Fourier transform of the power spectrum density according to the Wiener- Khinchin theorem, so the 
direct way is to derive the autocorrelation function. But this can cause problems when the data are 
irregularly spaced (Spangler et al. 1989). Regardless of the data sampling, it is easy to obtain the 
structure function, however. So we will take the structure function approach below. 

Given a sample of KM data, the logarithm of angular separations of RM pairs arc set to bins 
with equal lags. In each bin the differences of RM pairs are squared and then averaged to obtain the 
structure function. The standard deviations are taken as the uncertainty of each function value. 

Around each source we can derive the local RM average and dispersion. If a RM value is different 
from the average by three times larger than the dispersion, it is defined as "anomalous" . We will 
illustrate by simulation that such an anomalous RM can significantly affect the structure function. 

Let the RMs distribute in the region 0° < .t < 2° and 0° < y < 2° uniformly with grid intervals 
of 0.1°. Two samples of random RM values following Gaussian distributions with an average of zero 
and a dispersion a of 30 and 50 rad m~^, respectively, are generated with the Monte-Carlo method. 
Note that the average value does not affect the results. Theoretically the structure function should be 
flat with an amplitude 2a^ and an index 0. The structure functions from simulations are displayed in 
Fig. 1. We then fitted the structure functions with £>rm(^^) = A59°' and found that the results are 
consistent with theoretical expectations. 

As is obvious in Fig. 2, an anomalous RM value (1000 rad m"^^ indicated by a cross) can distort 
the real structure function to a different extent. With regard to actual RM data, such a RM can be 
either due to an enhanced local medium or it is intrinsic to the extragalactic source. One should pick 
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Fig. 1. Simulated RM distributions are shown in the lower panels, with filled and open circles representing positive 
and negative values, respectively. The sizes of symbols arc proportional to the square-root of RM values. The structure 
functions are plotted in the upper panels with the number of pairs in each bin marked. 



Table 1. Fitting parameters of the structure functions in the CGPS regions and the NGP region. 



Region 


Range 


(RM) 


CRM 


A 


a 


CGPSl 


82?42 < I < 96?23 


-237 


306 


118630±6901 


0.28 ±0.03 


CGPS2 


115?01 < I < 129?99 


-147 


117 


23016±1241 


0.10 ±0.03 


CGPS3 


130?12 < I < 146?51 


-106 


94 


17756± 958 


-0.003±0.027 


North Galactic Pole 


(6 > 70°) 


-3 


15 


548± 215 


-0.06 ±0.12 



up such RMs for further identification. Below we will only show the results of real observed RM data 
after discarding these anomalous RMs, as identifications have not yet been made. 

Wc took the RMs of the cxtragalactic sources (EGSs) in the CGPS region but excluded sources 
containing separate or unresolved components (marked by a, b and c in the catalog of Brown et al. 
2003). The data are divided into three longitude intervals: 82° < Z < 96° (CGPSl), 115° < I < 130° 
(CGPS2) and 130° < / < 146° (CGPS3). Wc also collected the RMs of EGSs in the North Galactic 
Pole (NGP) region with latitudes larger than 70°. Then we calculated the structure functions for these 
four regions, as shown in Figs. 3, 4, 5, and 6. The fitting results using DuuiSO) = ASO" are listed in 
Table 1, where Column 1 refers to the region. Columns 2 and 3 to the average and dispersion of the 
RMs. The amplitude and power index of each power law structure function are listed in Columns 4 
and 5. 



4 Discussion 



The electron density fluctuations are better known than the magnetic field fluctuations. Armstrong 
et al. (1995) have shown that the electron density fluctuation spectrum is Kolmogorov-like, but the 
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Fig. 2. Same as Fig. 1 for <t = 50 rad m ^ but with an anomalous RM of 1000 rad m ^ in the corner and in the 

center as indicated by a cross. 
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Fig. 3. The RM distribution in the CGPSl region (lower 
panel) and the calculated structure function (upper panel). 
The symbols are plotted as in Fig. 1, but anomalous RMs 
are plotted as hatched circles when positive or dotted circles 
when negative. 



Fig. 4. Same as Fig. 3 but for the CGPS2 region. 



Structure Function Studies for Turbulent Interstellar Medium 



29 



i5 - 



■A= 
.a= 


1 

17756±958 
-0.003±0.027 




« 


1 1 1 1 1 




t + 


» 

881 


2260 


3846 2099- 




74 








' U 












Totol: '(44 Discorded: § 



o -100 rod m ^ 



• oo, oOO <->^ ■ O " 
O o 




I (•) 



Fig. 5. Same as Fig. 3 but for the CGPS3 region. 
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Fig. 6. Same as Fig. 3 but for the NGP region. 



magnetic field fluctuation remains mysterious and totally unclear as mentioned above. We know that 
many physical processes in the ISM aS'ect the electron density and the magnetic field at the same time, 
such as supernova explosions. The magnetic field is often assumed to be frozen into the interstellar 
medium. So we can assume that there is no distinct discrepancy between the fluctuations of the electron 
density and the magnetic field. In this case the structure function can be written as £>rm ('5^') = 
^{66) + CSO"-'^ where C is a constant. 

In the three sub-CGPS regions, the average rotation measures (RM) = —237, —147 and —106 rad 
m~^, respectively, corresponding to the geometrical term contributions ^rm{S0) = IISO"^, 659"^ and 
366"^, which are definitely negligible when compared to the obtained amplitudes (A) of the structure 
functions in these regions. In fact, the indices from our results are much smaller than 2, which also 
means that the geometrical term plays a minor role in the structure function. Therefore the structure 
function can be simplified as Dy^m{SO) ~ C50°'~'^. 

The structure function Z?fim(<50) ~ CS6^^^ holds for 3D Kolmogorov fluctuations and -Dr,m('5^) ~ 
CSO"^/^ for 2D Kolmogorov fluctuations (Minter & Spangler 1996). It is clear that the Kolmogorov 
spectrum cannot account for our results. Actually, the indices of the structure functions in the three 
sub-CGPS regions are very small, close to zero. Let us try to discuss the nearly flat structure functions 
we got. 

From the results (Table 1) it is evident that the amplitudes of the structure function A decrease 
from CGPSl to CGPS3. This is probably caused by the extent of the ISM in our Galaxy. Prom our 
simulation, we know that A ~ 2a'j^j^j, where a^M is the dispersion of RM along the line of sight. 
Assuming that there are many fluctuation cells with typical scale I and a RM dispersion ui along the 
line of sight, the total RM dispersion can be represented by (t^^ = J^l'i where L is the path length 
to the edge of our Galaxy. Because L decreases from CGPSl to CGPS3, it is understandable that the 
fluctuation intensity gets smaller. 

Due to the disk structure of the Galactic electron distribution and the Sun's location at 8 kpc 
distance from the Galactic center, the RMs of EGSs near the pole region are little affected by the 
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Galactic interstellar medium, so that the structure function for the EGSs in the pole region should 
be very flat. Our result in Fig. 6 and the fitting parameters in Table 1 are consistent with Simonetti 
et al. (1984) and strengthen the fact that the RMs of EGSs in the pole region are very random, 
indeed mainly of intrinsic origin. Wc can also see that the fluctuation intensity in the NGP region is 
much smaller than that in the CGPS region, indicating a trend of a latitude-dependence of the ISM 
turbulence. 

5 Conclusion 

We have used the structure function to study ISM fluctuations. Simulation shows that an anomalous 

RM can significantly infiuence the structure function, so wc suggest that these large RMs should 
be carefully checked before performing a structure function analysis, or otherwise should be used 
with caution. The structure functions of the RMs of EGSs in the CGPS region have been calculated. 
The structure functions at SO > 1° can be fitted by a power law. The power indices are all nearly 
zero which cannot be interpreted by Kolmogorov-like ISM fluctuations. The fluctuation intensity is 
longitude-dependent. The flat structure function in the NGP region shows that the RMs in this region 
are almost intrinsic to EGSs. 
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